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Abstract 

A novel method is proposed to measure the Purcell effect by observing the 
current through a semiconductor quantum dot embedded inside a microcavity. 
The stationary current is shown to be altered if one varies the cavity length. For 
the double-dot system, the stationary current is found to show the interference 
feature (superradiance) as the inter-dot distance is varied. The amplitude of 
oscillation can be increased by incorporating the system into a microcavity. 
Furthermore, the current is suppressed if the dot distance is small compared 
to the wavelength of the emitted photon. This photon trapping phenomenon 
generates the entangled state and may be used to control the emission of single 
photons at predetermined times. 

1 Introduction 

Since Dicke proposed the idea of superradiance T , coherent radiation phenomena 
for atomic systems was intensively investigated (2J EH EJ 03 Ej ■ What Dicke found 
was that when the gas is in a particular state with the half number of molecules 
excited, the spontaneous emission rate of the whole gas is proportional to the square 
of the molecular concentration, provided that the gas volume dimension is small 
compared to the emitted photon wavelength. One of the interests in superradiant 
study lies in its close connection with the physics of laser emission. In some aspect, 
the superradiant phenomena appears somewhat simpler since one can neglect the 
pumping and relaxation mechanisms which are important in laser operation and 
consider only the evolution of atoms exclusively coupled to their own radiation 
field. 

In solid state physics, the exciton-polariton state is one of the limiting cases of su- 
perradiance. When a Frenkel exciton couples to the radiation field in a small system 
which contains iV lattice points, it represents one excited atom and N — 1 unexcited 
atoms in the others. According to Dicke's theory, the decay rate of the system will be 
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enhanced by a factor N. But as it was well known in a 3-D bulk crystal[7], the exci- 
ton will couple with photon to form polariton-the eigenstate of the combined system 
consisting of the crystal and the radiation field which does not decay radiatively. 
What makes the exciton trapped in the bulk crystal is the conservation of crystal 
momentum. If one considers a linear chain or a thin film, the radiative recombina- 
tion of excitons is fundamentally different from the three-dimensional case. In bulk 
crystal, radiative decay requires phonons or other translational-invariance breaking 
entities such as defects, impurities, or interfaces. Thin films or linear chains, how- 
ever, inherently break translational invariance. As a consequence, a radiative decay 
channel opens up, and the decay rate of the exciton is enhanced by a factor of X/d 
in a linear chain [H] and (X/d) 2 for 2D exciton-polariton|9l I10j. where A is the wave 
length of emitted photon and d is the lattice constant of the linear chain or the thin 
film. 

With the advances of microfabrication technologies such as molecular beam epi- 
taxy, it has become possible to fabricate various structures of microcrystals with 
fine quality and novel properties, such as quantum well, superlattice, quantum dot, 
quantum wire, and quantum ring. The exciton in a quantum well can exhibit the 
behavior between purely three dimensions and two dimensions. Many investiga- 
tions on the radiative linewidth of excitons in quantum wells have been performed 
[111 1121 1131 1141 1151 1161 117j . For lower dimensional systems, first observation of su- 
perradiant short lifetimes of excitons was performed by Ya. Aaviksoo et al. JE] on 
surface states of the anthracene crystal. Superradiance has also been discussed for 
one-dimensional (ID) Frenkel excitons in disordered aggregates. |19|l2lT] A. L. Ivanov 
and H. Haug^ST] predicted the existence of exciton crystal, which favors coherent 
emission in the form of superradiance, in quantum wires. Y. Manabe et al. 1221 
considered the superradiance of interacting Frenkel excitons in a linear chain. For 

o 

quantum dots of CuCl with radii R between 18 and 77 A, the superradiance of 
excitons was also observed by Nakamura et al. |23| . The decay rate was shown to 
be proportional to R 21 which confirms the theoretical prediction by Hanamura. |13j 
Similar works were also obtained by Itch et al. 24]. By using numerical simulation, 
Spano et al. |24| have also showed the effective coherent size (which is inversely 
proportional to the decay time in their definition) decreases as the temperature in- 
creases. Recently, superradiance of polaritons from a dimer to finite one-dimensional 
crystal has also been discussed by Dubovskii.[2l)] 

In reality, superradiance is accompanied by frequency shift, as pointed by Lee 
et aZ.[2Z|- However, the coherent frequency shift of an exciton has received fewer 
attention. One of the reasons is the difficulty arose from the divergent nature of 
frequency shift, both infrared and ultraviolet. Lee et al. [2Sj have solved the problem 
by using the method of renormalization for a system of two atoms, and applied to the 
case of excitons in a thin semiconductor film. [22j Recently, we have generalized their 
results to the quantum well systems. j29] The crossover behavior from 2D film to 3D 
crystal was also examined. It was found that both the decay rate and renormalized 
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frequency shift show oscillatory dependence on layer thickness. 

On the other hand, Purcell 30 predicted that the spontaneous emission rate, and 
thus the relaxation lifetime of an excited atom, would be altered if the atom was put 
in a cavity with dimensions comparable to the transition wavelength of the atom. 
The reason for modified spontaneous emission rate is that spontaneous emission can 
be viewed as stimulated emission, stimulated by vacuum-field fluctuations. Hence, 
the lifetime can be altered by modifying the photon density of states. 

In the last few years, experiments which have confined photons in low- 
dimensional semiconductor quantum structures by using optical microcavities were 
in progress. Researche on combined quantum confined 2D carriers and 2D photon 
states has been performed by using quantum wells embedded in planar microcavi- 
ties, yielding interesting physics both in the weak and strong carrier-photon coupling 
regimes. |31j Results from 2D carriers combined with ID (0D) photon states have 
also been studied |32| I33j . Introducing further degrees of carrier confinement in such 
microcavity structures is a natural trend in this field. Recently, the incorporation of 
quantum dots in planar and pillar microcavities have been reported. [31] A coupling 
between ID electron state and 2D photon states were also obtained by inserting an 
array of quantum wires into a planar microcavity. [35 Such systems lead to dramatic 
changes in the exciton-polariton dispersion 36] and exciton lifetime. [571 13*5] 

In this manuscript, an alternative way is proposed to observe the inhibited or en- 
hanced spontaneous emission by embedding a quantum dot inside a microcavity 2 1 . 
After the injection of an electron and hole into a quantum dot, a photon is generated 
by the recombination of the exciton. This process allows one to determine Purcell 
effect by measuring the current through the quantum dot. Similarly, by embedding 
two quantum dots inside the cavity and controlling the gate voltage of one of the 
dots, one can not only determine the superradiant effect by measuring the station- 
ary current, but also induces the entangled states, which is one of the fundamental 
requirements for quantum information processing . 

This paper is organized as follows. A brief description of superradiant and 
Purcell effects is reviewed in section II. In section III, we discussed the transport 
properties of a two level quantum dot embedded in a planar microcavity. Electrical 
measurements of superradiance and the generated entanglement are shown in section 
IV. Finally, overall conclusions are presented in the last section. 

2 A brief review of Superradiant and Purcell effect 
2.1 Spontaneous emission of two coupled atoms 

Spontaneous emission is one of the fundamental concepts in quantum mechanics 
that can be traced back to the early works of Albert Einstein. [22] In free space a 
two level atom interacts with a continuum of radiation field modes, which leads 
to an irreversible exponential decay of the excitation. In this subsection, the 
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Weisskopf-Wigner theory of spontaneous emission for two level atoms in free space 
is reviewed. |40| 

The total Hamiltonian for a single two level atom in the rotating wave approxi- 
mation (RWA) is given by 

h = + y, w&i&q + \) - E hD ^- + b ^+)> 

q q 

where a z = ||) (T| — U) (II and a + = ||) (|| , <r_ = | j) (|| are the Pauli matrices in the 
2x2 space, fiuj is the level spacing between the two levels, and L> q = ( 2e ^y ) 1//2 "e > -7? 

is the dipole coupling matrix element. Furthermore, u; q = c|q| , and 6 q creates a 
photon with wave vector q. 

In the interaction picture, the combined atom-field system is represented by 

|^)) = /o(t)|T;0) + ^/ G;q (t)|I;l g ), (2.2) 
q 

where |J.; l q ) represents the state in which the atom is in the ground state and the 
field mode q has one photon. 

Substituting Eq. (2.1) into Schrodinger equation and projecting into each state, 
we obtain 

/o(i) = exp[-(-r + *Aw)*], (2.3) 

where 



r = J ^ysiu - Wq ) 

(2.4) 



3TTEohc 3 



and 



1 P uj^ uP^ 

Aco = — £o hc 3 V / dwq—^—. (2.5) 

D7T J LJ — Cdq 

The spontaneous decay rate T gives the Einstein's A coefficient, and the frequency 
shift Alo represents the Lamb shift. [IT] 

In a system of two identical atoms interacting via common radiation field, the 
decay splits into a sub- and a superradiant channel. The Hamiltonian for two atoms 
interacting with the electromagnetic field reads 
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H — Hq + H p h + H ep h 

tt 1 /A A , 

H = -hu(a 1:Z + a 2 , z ) 
q 

H eph = ^>q(6^a ji _ + 6 q e-^ r ^ J -+), (2.6) 

q 

j=l,2 

where o^+and CTj i2 are the Pauli matrices in the 2x2 space of the upper/lower 
level It)- 7 , H} 3 of atom j, Huj is the level spacing between the two level, and £> q = 
( 2e ^y ) 1//2 "e > • Jl is the dipole coupling matrix element. Furthermore, u; q = c|q| , 

and 6 q creates a photon with wave vector q. 
One can define the so-called Dicke state, 



|Sb> = ^(ITD-IIT)) 

\Ti) = ITT) 

|T ) = -L(|U) + ||T)) 

\T-i) = ITT) • (2.7) 
Using this basis, one can easily calculate the matrix elements 



(Ti| <r jt± |Ti) = (Ti| aj, ± |T_i) = 0, j = 1, 2 
(Txl^+lTo) = (T |a j , + |r_ 1 ) = -^, j = l,2 

(To|a j)± |5o) = 0, j = 1,2 

(Ti|cti j+ |S ) = --j=, (Ti\<t 2 ,+ \S ) = -j= 

(Solk.+ lT-i) = -J=, <5o|a 2)+ |T_ 1 ) = --^. (2.8) 

This means that there are two transition rate T± for spontaneous of photons 
into a vacuum state, 

r ± («) = 2vr £ K± / q ' % - » q ), q = ~ c , (2.9) 
q 
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where we have defined a q = -D q e iqri and /3 q = D q e jqr2 . Evaluation of this expres- 
sion yields 

r ±( ,) = r[i ± ^], (2.10) 

where r oc q 3 is the decay rate of an isolated atom. Here, d = |ri — r2 1 is the 
distance between the two atoms. The appearance of two decay channels has been 
discovered by Dicke and observed by DeVoe and Brewer jl^l in a laser-trapped 
two-ion system. 

The time-dependence of the collective decay of two radiators is different from 
the decay of two single radiators. If we denote the occupation probabilities of the 
four levels by To(t), T_i(i), and So(t), the time dependence occupations is 

then governed by two decay rates r + and T_ : 



Tx = -(r + + r_)T! 
So = r_(Ti-sb) 

T = r+CTx-To) 

T_i = r + To + r_5o. 

The above equation can be solved easily, 



(2.11) 



Ti(t) = 




+r_)t 




So(t) = 


[e- T - 


-t_ e -(T. 


+r+)t] r _ 




r+ 




To(t) = 


[e- ^ ^ 


-t _ e -(r_ 


+ r +)*]r+ 




r_ 




T-i(t) = 


r_r_ 




fr+)*[(_i 



r_r. 



e r -*)ri 



(2.12) 



where initial conditions Ti(0) = 1, 5o(0) = To(0) = T_i(0) =0 have been assumed. 
If we consider the special case where T_ = and r + = 2T, this would correspond 
to the case qd — > 0, i.e. the wavelength of the emitted photon is much larger than 
the distance between the two radiators. Then, Eq. (2.12) reduced to 



T (t) = r+te- r +* 

So(t) = 

T_i(t) = 1 - e - r +*(l + r+t). (2.13) 
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The total coherent emission rate hit) at time t is the sum of the emission rates from 
Ti and T : 

hit) = £ r +e - r +*(i + r+t), r + = 2r, (2.14) 

where E$ is a constant with dimension energy. This is different to the incoherent 
sum 2ii(f) of the emission rates hit) from two independent radiators, which would 
give 

2hit) = 2E Te- rt . (2.15) 

2.2 Effect of cavity on the radiative decay of excitons in low di- 
mensional systems 

As we mentioned above, the electron-hole pair is naturally a candidate for examining 
spontaneous emission and Purcell effect. Let us first consider a Wannier exciton in 
a quantum ring with radius p ~ Nd/2ir, where d is the lattice spacing and N is 
the number of the lattice points. In our model, the circular ring is joined by the 
N lattice points, and we also assume the effective mass approximation is valid in 
the circumference direction. The state of the Wannier exciton can be specified as 
\u,n,m), where v is the exciton wave number, n and m are quantum numbers for 
internal structure of the exciton, and will be specified later. Here, v takes the value 
of an integer. The matter Hamiltonian can be written as 

Hex = ^ ^ E unm c\, nm Ci, nm , (2.16) 

vnm 

where ctnm and c unm are the creation and destruction operators of the exciton, 
respectively. The Hamiltonian of free photon is 

H ph = n <l' 2 + *?) V % K x h <iK>» ( 2 - 17 ) 

q>k> z 

where , k ,^ and & q '^A are the creation and destruction operators of the photon, 
respectively. The wave vector k' of the photon is separated into two parts: k' z is the 
perpendicular component of k' on the ring plane such that k' 2 = q' 2 + k' 2 . 
The interaction between the exciton and the photon can be expressed as 

H ' = D 4K™™ h K<i! c ™m + h - c -> ( 2 - 18 ) 

k' z nm q' 

where 



£>q'fci«mm = H^^p^J ~ 2 (2.19) 
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with eq'jf being the polarization of the photon and H, 
Eq. (2.19), 



(i) 



is the Hankel function. In 



x exp(iv(tp e< f* ))(-ih^-)w v (ip). 



d_ 

dtp' 



(2.20) 



is the effective transition dipole matrix element and F nm {ip e ) is the hydrogenic 
wavefunction in the ring. Here, m* e and m* h are, respectively, the effective masses 
of the electron and hole. 

The decay rate of the exciton can be expressed as 



"^f vnn 



2vr |Af 



k' unm 



(2.21) 



where C0q'k' z unm = E vnm /fi — cy q' 2 + k' 2 . The Wannier exciton decay rate in the 
optical region can be calculated straightforwardly and is given by 



e 2 h p 
m z c d 



where 



J2 



\eq>k' z \ ■ Xvnm\ dk' z dq , (2.22) 



/d 
d(pw*(<p - <p e )(-ih—)w v ((p). (2.23) 



From Eq. (2.22), one observes that the decay rate 7 vnm is proportional to p/d. 
This is just the superradiance factor coming from the coherent contributions of 
atoms within half a wavelength or so. In Fig. 1 we have numerically calculated the 
superradiant decay rate in v = 0, n = 0, and m = mode. In plotting the figure, 

o o 

we have assumed E unm /h = ko = 2ir/X, A = 8000vl, d = 5^4, 70 is the decay rate 
of an isolated atom, and for large radius, F* m is independent of p. The decay rate 
increases linearly with the increasing of radius when the radius is small. This linear 
regime agrees with Dicke's prediction: For one excited atom and N — 1 unexcited 
atoms in a small volume, the decay rate is enhanced by the factor of 1 x N. For 
large radius, the decay rate can approach ID limit (= 2^70) correctly. 

Quite recently, R. A. Romer and M. E. Raikh studied theoretically the exciton 
absorption shredded by a magnetic flux $ jl^. From their results, effects of magnetic 
flux on exciton wavefunction F* m can not be neglected in small radius limit, and 
may be examined from the variations of the decay rate. In the inset of Fig. 1, three 
curves of different flux $ are presented as functions of radius p. The dashed, solid, 
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Figure 1: Decay rate of a quantum ring exciton in large radius limit, i.e. F* m 
are assumed to be independent of p. The vertical unit and horizontal units are 
2§^7o and Ao, respectively. Inset : Effect of Aharonov-Bohm on the radiative decay 

of quantum ring exciton. The dashed( ), solid, and dashed-dotted(- •) curves 

correspond to $ = 0&o, 0.253>o> an d 0.5<l?o> respectively. In small radius limit, 
F* m depends strongly on radius p, and its influence on the decay rate is evident. 
The vertical and horizontal units here are 2^70 an d rrn g radius (in units of ao), 
respectively. 

and dotted curves represent the cases of <3? = O^o, 0.25<l?o> and 0.5<I>o> respectively. 
For <I> = O.S^O; the decay rate decreases as the ring radius becomes small but reaches 
the minimum point as p is about 0.25ao (where ao is the effective Bohr radius we 
assumed in ID limit). This is because the probability, for electron and hole to 
meet each other on the opposite side of the ring, increases with the decreasing of 
ring radius, while the coherent effect decreases with the decreasing of the radius. 
Therefore, there is a competition between these two effects as one decreases the 
radius. In Fig. 2, relative decay rate [7„„ m ($) — 7 I/ , im ( < J > = 0)] is plotted as a function 
of magnetic flux with different radius. The solid and dashed lines represent the 
cases of p = 1 ao and p = 0.5 ao , respectively. As expected, the larger the radius, 
the smaller the AB oscillation amplitude. Besides, the superradiant decay rate is 
most enhanced for <3? = 0.5<&0) an d the oscillation period is equal to <3?o = hc/e. 

We now consider a Wannier exciton in a quantum ring embedded in perfectly 
reflecting mirrors with cavity length L c . The decay rate of the quantum ring exciton 
can be expressed as 

|e q %-X,| 2 - (2-24) 

The numerical calculations of Eq. (2.24) are shown in the left panel of Fig. 3. As 
can be seen, the decay rate of a quantum ring exciton shows the enhanced peaks as 
the cavity length L c is equal to multiple half-wavelengths of the emitted photon. 




2~n 



P 

m?c 2 L r . d, 



i?( 1 )( v /(2vr/A)2-(vrn c /L c )2 /0 ) 
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Figure 2: Dependence of relative decay rate [y vnm {$) — 7„ nm ($ = 0)] on the mag- 
netic flux. The dashed and solid curves correspond to p = 0.5 ao and p = 1 ag, 
respectively. The vertical and horizontal units are 2^70 and universal flux quantum 
$o = hc/e, respectively. 



quantum ring quantum dot 
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Figure 3: Left panel : Decay rate of a quantum ring exciton in a planar micro- 
cavity with radius p = \/2ir. The horizontal and vertical units are (A/2) and 
^Srrfic^d l e q' fc z ' Xv\ , respectively. Right Panel : Similar case for a quantum dot 
exciton. 



However, if one considers a quantum dot exciton inside the microcavity, the 
decay rate reads 

7 « £ -Sr^(( 27r / A ) 2 " (™c/L c f) ■ xf , (2-25) 

TYt C J-j f 

where 9 is the step function. The numerical calculations are presented in the right 
panel of Fig. 3. One can see from the figure, there is no enhanced peak with the 
increasing of the cavity length. This is because the angular momentum (translational 
momentum) of the exciton in a quantum ring is conserved in circular direction, while 
the crystal symmetry is totally broken in a quantum dot. Due to the modification 
of the density of states of the photon in the microcavity, the decay rate of the 
exciton shows enhanced peaks in ID systems [38] and zigzag structure in 0D quantum 
dot. One also notes that such kind of peak maybe a useful feature to realize the 
Aharonov-Bohm effect for an exciton in a quantum ring. Generally speaking, the 
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excitonic AB oscillation is very small and hard to be measured. However, if one can 
incorporate the quantum ring inside the planar microcavities, the AB oscillation 
may be enhanced at these peaks. 

3 Current through one quantum dot and Purcell effect 

We now consider a quantum dot embedded in a p-i-n junction which is similar to 
the device proposed by O. Benson et a/|44|. The energy-band diagram is shown in 
Fig. 4. 



p-GaAs 


InAs 






QD 




Electron subband 






exciton 








Electron 



^^^^^^ 



hole 
■ tunneling 




Figure 4: Energy-band diagram of the p-i-n junction. 



Both the hole and electron reservoirs are assumed to be in thermal equilibrium. 
For the physical phenomena we are interested in, the fermi level of the p(n)-side 
hole (electron) is slightly lower (higher) than the hole (electron) subband in the dot. 
After a hole is injected into the hole subband in the quantum dot, the n-side electron 
can tunnel into the exciton level because of the Coulomb interaction between the 
electron and hole. Thus, we may assume three dot states 



|0> = 10,/!) 
\U) = \e,h) 

\D) = |0,0) (3.1) 

, where |0, h) means there is one hole in the quantum dot, \e, h) is the exciton state, 
and |0,0) represents the ground state with no hole and electron in the quantum 
dot. One might argue that one can not neglect the state |e, 0) for real device since 
the tunable variable is the applied voltage. This can be resolved by fabricating a 
thicker barrier on the electron side so that the probability for an electron to tunnel 
in advance is very small. Moreover, the charged exciton and biexcitons states are 
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also neglected in our calculations. This means a low injection limit is required in the 
experiment |45|. We can now define the dot-operators ny = \U) (U\ , ud = \D) {D\ , 
p = \U) (D\ , sir = |0) (U\ , so = |0) (D\. The total hamiltonian H of the system 
consists of three parts: the dot hamiltonian, the photon bath, and the electron 
(hole) reservoirs: 

H = Hq + Ht + Hy 

Hq = e\jn\j + + Hp + H res 

H T = J2g(D k b)£ + Dtb$) = g$X + tfxi) 

k 

Hp = QJkbjbk 

k 

H v = ^{V^ih + W^d\s A D + c.c.) 
q 

Hres = ^Jctcq + ^ejdtdq. (3.2) 

q q 

In above equations, is the photon operator, gD^ is the dipole coupling strength, 
X = Yl,k Dkb\ , and c q and d q denote the electron operators in the left ad right 
reservoirs, respectively. Here, g is a constant with a unit of the tunneling rate. 
The couplings to the electron and hole reservoirs are given by the standard tunnel 
hamiltonian Hy, where Vq and couple the channels q of the electron and the 
hole reservoirs. If the couplings to the electron and the hole reservoirs are weak, 
then it is reasonable to assume that the standard Born-Markov approximation with 
respect to these couplings is valid. In this case, we will derive a master equation 
from the exact time-evolution of the system. 

A 

In interaction picture, time evolutions of arbitrary operators O and X are defined 

by 

5(t) = e iHot Oe- lHot , X t = e iHot Xe- iHot . (3.3) 
Furthermore, for the total density matrix H(t) which obeys the Liouville equation 

E(t) = e- iHt E t=0 e iHt , (3.4) 

and we also define 

S(t) = e iHot E(t)e- iHot . (3.5) 

A 

The expectation value of any operator O is given by 

(0) t = Tr(~(t)0) = Tr(E(t)0(t)). (3.6) 
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We therefore have 



nu(i) = nu, no{t) = no 
p(t) = pe i£t X u p^(t)=pe-^xt 

£ = £{/- ££>■ (3.7) 

The equation of motion for E(t) becomes 

i^E(t) = [H T (t) + Hy(t),E(t)]. (3.8) 

This can be written as 



j t m = -i[H T ( t )M)]-i[H^(t),m] 

- -i[H T (t),E(t)]-i[Hy(t),E ] 

- f\t'[Hy(t),[H T (t')+7^(t'),E(t')}}. (3.9) 
J o 

Now, we define the effective density operator of the dot plus photons, 

pit) = Tr res E(t) (3.10) 

as the trace of E(t) over electron reservoirs. The trace Tr res over the terms linear 
which are in Hy vanishes, therefore, 

j t p(t) = -i[H T (t),p(t)} -Tr res j\t\Hv(t),[Hy(t'),E(t')]}. (3.11) 

As can be seen from the above equation, the last term is already second order in 
Hy, we can approximate 

S(0 « R p(t'), (3.12) 

where Rq is the equilibrium density matrix for the two electron reservoirs. Working 
out the commutators and using the time evolution of the electron reservoir operators 

c~ (t) = d^(t) = e-< l d qi (3.13) 

the master equation becomes 



214 



Y.N. Chen, D.S. Chuu and T. Brandes 



p(t) = po-i f dt'[H T (t'),p{t')] 
Jo 

-T L f dt'{^(t')^(t')p(t') - 2^j\t')p(t')^j(t')} 
Jo 

-T L f dt'{p(t')^(t')^(t')} 
Jo 

-Tr f dt'{Zh\t')^(t')p{t')} 

Jo 

-T R f dt'{-2^(t')p(t')^\t')+p(t')^\t')^(t')}, (3.14) 
Jo 

where T L = vr £ q V^S{e v - e£) and T R = tt £ q W*6(e D - e*). 

Multiplying Eq. (3.14) by njj, no, p, and p , respectively and performing the 
trace with the three dot states in Eq. (3.1), one obtains 



(nu) t - (nu) = -ig I dt'{(p) t , - (p t ) } + 2T V I dt'(l - (nu) t , - (n D ) t , 



{n D ) t ~ (n D } = -igJ Q dt'{{p) t , - (p^J - 2T D J dt'(n D ) 



(p) t -(p)t = -r D J\t'e^-^(x t x},p(t')) tl 

-ig J* dt'e^HlnuXtXl)^ - (n D x\,X t ) t ) 



(pt) t - { pf t = -r D jT dt'e-^-v (pt(t>)x t ,xl) t/ 

+ig J^dt'e'^-^iinuX^xl)^ - (n D x\x t ,) t ), (3.15) 

where e = etj — Ed is the energy gap of the quantum dot exciton. Here, p(t') = 
pe t6t X t r, and X t > denotes the time evolution of X with H p . The expectation value 

A 

(p^) t describes the decay of an initial polarization of the system and plays no role 
for the stationary current. Therefore, we shall assume the initial expectation value 

of p vanishes at time t = 0. 
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As can be seen from Eqs. (3.15), there are terms like yiijX t x\^j ^ which con- 
tain products of dot operators and photon operators. If we are interested in small 
coupling parameters here, a decoupling of the reduced density matrix pit') can be 
written as 

pit') « p° ph Tr ph p(t')- (3-16) 
By using above equation, we obtain 

Trip^nuXtX],) « {nu) t , (*t*J) (3-17) 

and correspondingly the other products of operators can be obtained also. For spon- 
taneous emission, the photon bath is assumed to be in equilibrium. The expectation 

value (x t x},^ = C(t-t') is a function of the time interval only. We can now define 

the Laplace transformation for real z, 



/•oo 

C e (z) = / dte~ zt e i£t C{t) 
Jo 



oo A 

nu(z) = I dte- zt ( nu ) t etc., z>0 (3.18) 



and transform the whole equations of motion into z-space, 

nu(z) = -i%(z)-p*(z))+ 2-^(1/ z- nu (z)-n D (z)) 

n D {z) = |(p(z) -p*(z)) - 2^j-n D (z) 
p(z) = -ig{nu(z)C £ (z) - n D (z)Ct £ (z)} -T D p(z)C £ (z) 

p*(z) = ig{n u (z)C* £ (z)-n D (z)C- £ (z)}-T D p*(z)C* £ (z). (3.19) 

These equations can then be solved algebraically. The tunnel current / can be 

defined as the change of the occupation of nu and is given by I = ig(p — p ), where 
we have set the electron charge e = 1 for convenience. The time dependence of the 

A 

expectation value (I) t can be obtained by solving Eqs. (3.19) and performing the 
inverse Laplace transformation. For time t — ► oo, the result is 

A _ 2g 2 T u T D B 

[ 9 2 T D B + [g 2 B + T D + 2~fTl + ( 7 2 + n 2 )T 3 D ] 

B = 7 +( 7 2 + ft 2 )r D , (3.20) 

where g 2 VL and g 2 j are the exciton frequency shift and decay rate, respectively. The 
derivation of the current equation is closely analogous to the spontaneous emission 
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of phonons in double dots |46| . in which the correlation functions (^XtXj,J is given 
by the electron-phonon interaction. 

As can be seen from Eq. (3.20), the stationary current through the quantum dot 
depends strongly on the decay rate 7. The tunnel currents of a quantum dot inside 
a planar microcavity is numerically displayed in Fig. 5. In plotting the figure, the 
current is in terms of 100 pA, and the cavity length is in units of Ao/2, where Ao 
is the wavelength of the emitted photon. Furthermore, the tunneling rates, Tjj and 
Td, are assumed to be equal to O.270 and 70, respectively. Here, a value of 1/1. 3ns 
for the free-space quantum dot decay rate 70 is used in our calculations [37j. Also, 
the planar microcavity has a Lorentzian broadening at each resonant modes (with 
broadening width equals to 1% of each resonant mode) As the cavity length 

is less than half of the wavelength of the emitted photon, the stationary current is 
inhibited. This is because the energy of the photon generated by the quantum dot 
is less than the cut-off frequency of the planar microcavity. Moreover, the current 
is increased whenever the cavity length is equal to multiple half wavelength of the 
emitted photon. It represents as the cavity length exceeds some multiple wavelength, 
it opens up another decay channel abruptly for the quantum dot exciton, and turns 
out that the current is increased. With the increasing of cavity length, the stationary 
current becomes less affected by the cavity and gradually approaches to free space 
limit. 

0.3 

0.25 

< 

o 0.2 

o 

v 0.15 
c 

QJ 
M 

£ 0.1 

u 

0.05 



12 3 4 

Cavity Length [A /2] 

Figure 5: Stationary tunnel current, Eq. (3.20), as a function of cavity length L c . 

A 

The vertical and horizontal units are 100 pA and Ao, respectively. Inset : (J) as a 
function of exciton energy gap e. The cavity length is fixed to Ao/2. The current is 
in units of 100 pA, while the energy gap is terms of 2/ic/Ao- 

To understand the inhibited current thoroughly, we now fix the cavity length 
equal to Ao/2 and vary the exciton energy gap, while the planar microcavity is now 
assumed to be perfect. The vertical and horizontal units in the inset of Fig. 5 are 
100 pA and 2/ic/Ao, respectively. Here, Ao is the wavelength of the photon emitted 
by the quantum dot exciton in free space. Once again, we observe the suppressed 
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current as the exciton energy gap is tuned below the cut-off frequency. The plateau 
features in the inset of Fig. 5 also comes from the abruptly opened decay channels 
for the quantum dot exciton. From the experimental point of view, it is not possible 
to tune either the cavity length or the energy gap for such a wide range. A possible 
way is to vary the exciton gap around the first discontinuous point 2/ic/Ao- Since 
the discontinuities should smear out for the real microcavity, it is likely to have a 

A 

peak if one measures the differential conductance d{I)/de as a function of energy 



4 Current through the double-dot system and the in- 
duced entanglement 

Now, we consider two spatially separated quantum dots incorporated inside the p-i- 
n junction. The novel feature here is the dissipative creation of entanglement over 
relatively large distances, and its readout via the stationary current. The device 
structure is shown in Fig. 6. 



Figure 6: Proposed device structure. Two InAs quantum dots are embedded in 
a p-i-n junction. Above dot 2 is a metal gate, which control the energy gap and 
orientation of the dipole. 

One of the obstacles in measuring superradiance between the quantum dots 
comes from the random size of the dots which result in a random distribution of 
energy gap and thus diminishes the coherent radiation. This can be overcome by 
constructing a gate voltage over one of the quantum dots. The energy gap and the 
orientation of the dipole moments of one of the quantum dots, thus can be controlled 
well. 

After a hole is injected into the hole subband in the quantum dot, the n-side 
electron can tunnel into the exciton level because of the Coulomb interaction be- 
tween the electron and hole. In our calculation, we also neglect the Forster pro- 
cess which may have some influences on the results if the two dots are close to 
each other. The validity of this assumption will be discussed later. Thus, we 



gap e. 




p-ohmic 



dot 
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may assume four dot states |0) = |0, h;0, 0), \Ui) = |e, h;0, 0), 1 1/2) = \0,0;e,h), 
and |D) = |0, 0; 0, 0) , where |0,/i;0, 0) means there is one hole in dot 1 and 
|0, 0; 0, 0) represents the ground state with no hole and electron in the quantum 
dots. The exciton states |e, h;0, 0) (in dot 1) can be converted to \0,0;e,h) (in 
dot 2) through the exciton-photon interactions. By transforming \U\) and ^2) into 
Dicke states: \S ) = (|*7i) - \U 2 )) and \T ) = -±= (\Ui) + \U 2 )), we can now de- 
fine the dot-operators ns = \So) (Sol ,n% = \Tq) (Tq\ , no = \D) (D\ , p s = \So) (D\ , 
pr = |T ) (D\ ,su, = ^(|0) (S \ + |0) (T |), Sf) = |0) (D\. The exciton-photon cou- 
pling is described by an interaction Hamiltonian Ht- 

Ht = ^^Wk[P5(l + e lk - r ) 
k ^ 2 

+p T (l-e^ r )] + c.c.} 
= g(p A s X s +Ps r xJ+p A T Xr+p T t Xr t ), (4.1) 

where r is the position vector between two quantum dot, X$ = SkU + e * k r )^k^k 5 
and Xt = X^k(l ~~ e * k r )Ac^k- The dipole approximation is not used in our calcu- 
lation since we keep the full e lkr terms in the operators X$ and Xt- Following the 
derivations in previous section, one can also derive a master equation for this double 
dot system. The equations of motion can be expressed as 



A A ft A , A v 

(n a ) t -(n a ) Q = -W J o dt'{(p<r} t > ~ \Pl) t ,} 

+IV I dt'(l - (n s } t , - (n T ) t , - (n D ) t ,) 
Jo 

A A ft i A , A > A , A > 

(n D ) t ~ (nD} = -ig J Q dt'{{p s ) t , - (plj^ + {p T ) t , - [Pt/J 

ft A 

-2T D / dt'(n D ) t , 
Jo 

A A ft 

(ps) t ~ (Ps) t = 



\o 



(Pr) t - (pr) t = 



-r D ^ df e «(*-*') (x t 4^(i')) t , 

j\t'e^-^{^sX t x\) ti - (n D x\,X t ) t ) 
f dt'e^{(n A T X t xl) ti - (n A 4l ( ) f) }, (4.2) 



-W 1 
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where the index a = S or T. 

Similarly, the tunnel current / can be defined as the change of the occupation of 

A ^ A A t ^ 

no and is given by I = ig ^^(pa —pa )• The expectation value (I) t can be obtained 
in the limit of t —* oo and reads 

A 



7- + 7+[l + ^-(g 2 /^D +g 2 /Tu + T D ) 



(4.3) 



where <7 2 7+ and g 2 j- are the superradiant and subradiant decay rate of the exciton, 
respectively. 

The corresponding decay rate for superradiant and the subradiant channels is 
given by 



g 2 i± 



sin(2vr(i/A ) 
7o(! ± , /x ), 



(4.4) 



2-Kd/Xo 

where d is the inter-dot distance and 70 is the exciton decay rate in a quantum dot. 
To display the dependence of the stationary current through the quantum dot on the 
dot distance d, we present the results of two identical quantum dots in Fig. 7. As 



. 3 

< 0.25 

a 

o 

S 0.2 

£ 0.15 

u 
u 

P 0.1 

u 

.05 




0.5 1 1.5 2 2.5 3 

dot distance 



0.6 0.8 1 

dot distance 



Figure 7: Stationary tunnel current as a function of dot distance d. The interfer- 
ence effect is seen clearly (inset) by incorporating the system inside a rectangular 
microcavity. The vertical and horizontal units are 100 pA and Ao, respectively. 

shown in Fig. 7, the current is suppressed as the dot distance d is much smaller than 
the wavelength (Ao) of the emitted photon. This corresponds to the trapping state 
in the two-ion system. As long as we choose only one of the dots to be coupled to 
reservoirs, the generated photon is reabsorbed immediately by the other dot and vice 
versa. The current is then blocked by this exchange process. For small rates limit 
(<? 2 7±) one can approximate Eq. (4.3) by I w 4[l/g 2 7_ + l/g 2 ^/ + ]~ 1 . The rates Tjj,d 
drop out completely and the current is only determined by the (smaller) radiative 
decay rates. In this approximate form, the current looks identical to the expression 
for the conductance G cc [1/Tl, + l/T^] -1 from a left lead through a single level 
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to a right lead with tunnel rates This implies that the superradiant and the 

subradiant channel are in series (and not in parallel) in this limit. This is because 
once the exciton is formed in dot 1, time evolution of this state is proportional 
to e~ 9 7+4 + e~ 9 7 ~* not e~ 9 (t++T-)* ; 28^. It means the two decay channels in our 
system are not in parallel. For long time behavior t — > oo and 7+ >> 7_, the 
function e~ 9 7+i + e~ 9 7 ~* approaches the limit of e~ 9 7 ~*, which is identical to the 



same limit of the function e (in series). 

Similar to the two-ion superradiance 02]) the current also exhibits oscillatory 
behavior as a function of dot distance. To observe the interference effect clearly, one 
may incorporate the system inside a microcavity since semiconductor cavities with 
strong electron-photon coupling have been realized experimentally by, e.g., Gerard et 
al. 48 ]. Reduction of the allowed /c-state is expected to increase the magnitude of the 
oscillation. For example, if the system is placed inside a rectangular microcavity with 
length Ao, the decay rate for the two channels can be worked out straightforwardly: 



The stationary current is plotted in the inset of Fig. 7, where a perfect (lossless) 
cavity is assumed. As we mentioned above, the amplitude of oscillation is larger than 
that in free space. However, the oscillation period is not half of the wavelength, but 
A /\/2. This is because the interference term is only influenced by the wave vector 
in the unconfined direction. Excluding the contributions from fundamental cavity 
modes, the effective wave vector can be expressed as 



The oscillation period of the decay rate and the current is therefore increased by a 
factor of y/2. 

In Fig. 8, we plot the expectation value of ns (tit) as a function of the dot 
distance. The maximum entangled state (|So)) is reached as d « Ao- This is 
remarkable as the steady state is independent of the initial state. The entanglement 
is induced by the cooperative decoherence in the system. In a recent paper by 
Schnider et al.\4§\, the authors consider the behavior of an ion trap with all ions 
driven simultaneously and coupled collectively to a heat bath. They also found that 
the steady state of the ion trap can exhibit quantum entanglement. However, the 
concurrence of their system is below the value of unity (maximum entanglement). 
On the contrary, in our system the maximum entangled state can be generated 
by tuning the band gap of dot 2 (linear stark effect), i.e. control the on/off of 
the superradiance. Another advantage of our scheme is shown in the inset of Fig. 
8. If the double-dot system is incorporated inside a rectangular microcavity, the 
maximum entangled states repeat as a function of inter-dot distance. This means 
even for remote separation, the entanglement can still be achieved. The reason 





= ™ lie 

7T 



,i27rd/(v / 2Ao) 2 



(4.5) 




(4.6) 
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can be attributed to that the creation of entanglement in our model is governed 
by the interaction with a common heat bath [SO], while conventional creation of 
entanglement depends on the direct interaction between two subsystems |51j. When 
two dots are coupled to the common photon fields, the collective decay process drives 
the system into the entangled states. The novel feature of the effect predicted here 
is that entanglement in fact can be controlled electrically (without applying a laser 
field) and read out in the form of a transport property, i.e., the electron current (as 
a function of the dot distance or, alternatively, the cavity length). 




0.5 1 1.5 2 2.5 3 

dot distance 



Figure 8: Occupation probability of the entangled states ns (solid line) and ny 
(dashed line). The inset shows the results inside a rectangular microcavity. 

Another possible application of this effect is that by tuning the coherence of the 
dots, one can control the emission of single photon at predetermined times, which is 
important for the field of quantum information technology. One might argue that for 
small inter-dot distance the Forster process may play some role in our system 52 ; 
nevertheless, this only causes small energy splitting between state \So) and \Tq). 
Comparing to the large energy difference in the III-V semiconductor material, its 
effect on the decay rate g 2 ~/± is negligible. 

A few remarks about the problem of dissipation should be mentioned here. The 
coherence of the quantum states is a fundamental issue in quantum physics and 
decoherence caused by phonons or imperfections may destroy the unitary quantum 
evolution. In our proposals, decoherence due to interaction with other bosonic 
excitations (phonons and electron-hole pairs in the leads) is inevitable but can in 
principle be (partly) controlled by variation of the dot energies, or control of the 
mechanical degree of freedom [53]. In addition, scattering due to impurities are 
negligible since there is no interdot transport in our system. 

5 Conclusions 

In summary, we have proposed a new method of detecting superradiant and Purcell 
effects in semiconductor quantum dots. By incorporating a quantum dot between 
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a p-i-n junction, the Purcell effect on stationary tunnel current can be examined 
by the variations of cavity length or exciton gap. For the double-dot system, the 
superradiant effect on the stationary tunnel current can be examined by tuning the 
band gap of the quantum dot. The interference effects between two dots can be 
seen more clearly by incorporating the system inside a microcavity. The oscillation 
period of the decay rate and current is also increased because of the microcavity. 
Moreover, the maximum entangled state is induced as the inter-dot distance is much 
smaller than the wavelength of the emitted photon. Our model provides a new way 
to generate the entanglement in solid-state systems and maybe useful in future 
quantum information processing. 
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